Phonon-assisted tunneling and two-channel Kondo physics in molecular junctions. 
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The interplay between vibrational modes and Kondo physics is a fundamental aspect of transport 
properties of correlated molecular conductors. We present theoretical results for a single molecule 
in the Kondo regime connected to left and right metallic leads, creating the usual coupling to a 
conduction channel with left-right parity ("even"). A center-of-mass vibrational mode introduces 
an additional, phonon-assisted, tunneling through the antisymmetric ("odd") channel. A non- 
Fermi liquid fixed point, reminiscent of the two-channel Kondo effect, appears at a critical value of 
the phonon-mediated coupling strength. Our numerical renormalization-group calculations for this 
system reveal non-Fermi-liquid behavior at low temperatures over lines of critical points. Signatures 
of this strongly correlated state are prominent in the thermodynamic properties and in the linear 
conductance. 
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I. INTRODUCTION 

Ground-breaking experimental results in single molec- 
ular transistors during the last decade^ have greatly ex- 
panded the field of molecular electronics, opening sev- 
eral possibilities for technological applications and inves- 
tigations of fundamental aspects of the physics of these 
devices. By now it is established that strong correla- 
tion effects play a key role in the electronic transport 
through these systems, as evidenced by the observation 
of the Kondo effect 2 - in both break junctions 3 -^ and STM 
setups^ A clear understanding of the mechanisms in- 
volved in the emergence of the Kondo effect in molecular 
systems is thus of primary importance. 

One possibility to advance our knowledge of these de- 
vices is by establishing analogies with the well-known 
transport properties of semiconductor quantum dots in 
the Kondo regime. This, however, proves to be a chal- 
lenging approach for several reasons: molecule-leads cou- 
plings are very sensitive to the particular configurations, 
charging energies are significantly larger and, more im- 
portantly, deformations and vibrational modes in the 
molecule play an active role in transport £ 

This variety of competing effects also brings theoreti- 
cal challenges, such as the interesting prospect of investi- 
gating the interplay between vibronic states and Kondo 
physics. Different studies have investigated the effect 
of electron-phonon couplings in the charge degrees of 
freedom of the molecule, affecting the exchange corre- 
lations leading to the Kondo effect jLSiS This issue has 
also been highlighted in recent experiments reporting 
anomalous behavior in the Kondo transport which have 
been attributed to the "dressing" of the local energies by 
Holstein-like phononsi 7 - 

In addition to these local effects, considerable atten- 
tion has been given to effects of vibrational modes in 
the tunneling from the molecule to the leads. Phonon- 
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assisted couplings by "breathing 
mass" 10 ' 12 i 13 ' 14 molecular modes create additional corre- 



lations with the electrons in the leads. Such phonon- 
mediated tunneling processes will, in general, lead to 
novel features in the transport propertie s 10 ' 11 ! 12 ? 13 ! 14 and 
can be experimentally probed by conductance measure- 
ments. 

We address this subject in the present work by investi- 
gating a two-channel Kondo (2chK) effec t 15 ' 16 in molec- 
ular systems with "center-of-mass" vibrational modes. 
Two-channel Kondo physics, originally investigated in 
the context of heavy- fermion materials) 15 ' 17 has been an 
active topic in the area of nanostructures. In semicon- 
ductor quantum dots, several theoretical predictions^ 
and a recent experimental observation^ of the 2chK have 
highlighted the renewed interest in such strongly corre- 
lated states. In addition, phonon-assisted 2chK behavior 
has been predicted in effectively nonintcracting systems 
(e.g., metallic carbon nanotubes 2 ^ and metallic break 
junctions 2 ^) coupled to vibronic states with Kondo-likc 
correlations appearing in orbital (pseudospin) degrees of 
freedomi 20 1 21 

In this paper, we consider the Kondo regime of a singly 
charged molecular level connected to metallic leads, fully 
including both electron-electron and electron-phonon in- 
teractions, as well as phonon-assisted tunneling processes 
arising from a center-of-mass vibronic mode. Our nu- 
merical renormalization-group (NRG) calculations show 
that the presence of the extra, phonon-mediated, con- 
ductance channel leads to non-Fermi-liquid (NFL) be- 
havior at low temperatures, with prominent signatures 
in the thermodynamic properties and in the linear con- 
ductance. The 2chK fixed point occurs over lines of crit- 
ical points covering a wide range of parameters, both 
away from particle-hole symmetry and in the presence of 
deformation-induced charge-phonon couplings. 

The paper is organized as follows: the model is pre- 
sented in detail in Sec. |TT] and a discussion on the 2chK 
regime and the dependence of the critical parameters is 
given in Sec. IIII1 In Sec. IIV1 we discuss the NFL sig- 
natures in the transmission phase shift and in the con- 
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FIG. 1: (color online) Schematic illustration of the phonon- 
assisted tunneling process: the coupling of a molecular level to 
left and right metallic leads is modulated by a center-of-mass 
vibrational mode. 

ductance across the junction. We give our concluding 
remarks in Sec. [Vj 

II. MODEL 

We consider a molecular complex (e.g., an organo- 
metallic compound, Ceo, etc.) in the Coulomb blockade 
(CB) regime connected to metallic leads (for instance, in 
a metallic break junction setup). Kondo correlations ap- 
pear as the molecule is tuned into a CB valley with an 
odd number of electrons by a plunger gate voltage. We 
focus on the low-bias regime, in which electronic trans- 
port is dominated by a singly-occupied molecular level of 
energy e d (measured from the Fermi energy in the leads 
and tunable by the gate voltage) with a charging energy 
U arising from the electron-electron interactions within 
the molecule. 

In our model, the molecule is connected to left (L) and 
right (R) leads by tunneling couplings proportional to the 
overlap between the wavefunctions of the molecular level 
and the (s-like) metallic states in the leads. More impor- 
tantly, we consider the effect of phonon-mediated cou- 
pling through a center-of-mass vibrational mode, as illus- 
trated in Fig. [TJ For small displacements, the molecule- 
lead tunneling can be effectively written as V^^(l±ax) 
where x is a displacement operator in the direction of the 
motion and a is a system-specific parameter, being essen- 
tially proportional to the ratio between the overlap length 
and the oscillation amplitude of the vibrational mode. 
Additionally, deformations in the chemical bonding (also 
illustrated in Fig. [l} will, in general, introduce an effec- 
tive coupling of the center-of-mass oscillating mode and 
the charge state of the molecule. We model this by a 
Holstcin-type electron-phonon term with coupling A. 

The full Hamiltonian is given by 

H = Hm + H Loads + -^Dot-Leads , (1) 

where 

H M = £d"nda + Un d] n di + 



A(l — nd){a + a 1 ) + ujoa'a , 
-^Dot-Leads = ^ Vl (1 - ax) d\c Lka + h.c. 

k 

+V R (1 + ax) d^CBka + h.c., 

#Lcads = ^ ££fcC£ kcr C£kcr ■ (2) 
k,e=L,R 

In the above, (d a ) and cJ ko . (c^ko-) are fermionic 
operators that create (destroy) electrons with spin a in 
the molecule and leads, respectively, (n d a = d^d^ is the 
electron number operator), luq is the frequency of the 
local center-of-mass phonon mode, with a' (a) being the 
phonon operators (x = a+a^). We assume the wide band 
limit and k— independent dot-lead couplings (14k = Vt). 

Hamiltonian JT]) can be written as an Anderson im- 
purity model coupled to two independent fermionic 
channels. Defining symmetric ("even") and anti- 
symmetric ( "odd" ) combinations of the electronic 
operators in the left and right leads c e ( )k CT = 
{V R c Rka ± V L c Lka ) /2i/V^ + V|, the ^Dot-Leads term in 
([1]) becomes: 

Hoot -Leads — V rf^c e k CT + a (a + a^) d^c kcr + h.c. , 

k.tr 

(3) 

where V = 2\/V£ + VS. For a/0, a phonon-mediated 
coupling to the odd channel is present. As we shall 
see, this has important consequences in the physics of 
the ground-state of the system. Notice that the odd- 
channel coupling is present even for molecules not sym- 
metrically coupled to the leads (i.e., Vl ^ Vr), a 
more likely configuration in experiments. For a = 0, 
this term vanishes and Hamiltonian ([T]) corresponds to 
the single-channel Andcrson-Holstcin model, previously 
investigated by NRG^ and analytical renormalization- 
group methods £ 

The eigenstates of the Hamiltonian ([1]) can be labeled 
by total charge and SU(2) spin symmetry (Q,S). At 
the molecule site, states are also labeled by the number 
of phonons m ((a^) |0) = |m)) Notice that both the 
electron-phonon (cx A) and the phonon-assistcd tunneling 
(oc a) terms couple m) and |m±l) states. This last term 
couples only states with a difference of one electron in 
the odd channel and, thus, commutes with a generalized 
parity operator, defined by P = (— i) m +Q°+ 1 where Q a 
is the total charge in the odd channel. Therefore, for 
A = 0, the Hamiltonian |T]) has an additional 0(1) parity 
symmetry (which is lost for A ^ 0, as the electron-phonon 
term will couple states with different parity). 

We solve the two-channel problem with Wilson's 
numerical rcnormalization-group (NRG) technique^ 
adapted to include the phonon degrees of freedom 4 As 
it is standard in the NRG method, we map Eq. |T]) into a 
(two-band) tight-binding Hamiltonian by performing log- 
arithmic discretizations of the continuum spectra in the 
even and odd fermionic channels. In the calculations, 
we use a discretization parameter A = 3 and kept up to 



FIG. 2: (color online) Impurity contribution to the entropy 
'S'imp at the particle-hole symmetric point ed = —U/2 (a,c) 
and at Ed = — 0.9C/ (b,d). In both cases, an unstable 2chK 
fixed point for which S'imp — > k,B log 2 1 / 2 as T — > is reached 
at q = q c . Panels (a) and (b) show S'imp versus T/D for fixed 
a while (c) and (d) show Si mp versus a for fixed temperatures. 



2500 states in the NRG iterations, which proves to be 
adequate for the calculation of thermodynamic proper- 
ties. We have used a cutoff of N p ^ = 9 in the maximum 
number of phonons and checked for convergence of the re- 
sults with iVph. We assume a constant (metallic) density 
of states po in the leads with bandwidth D. 



III. TWO-CHANNEL KONDO PHYSICS 

From the NRG spectra, we calculated thermodynamic 
properties for this system. A particularly revealing quan- 
tity is the contribution to the total entropy coming from 
the "impurity" degrees of freedom (in the present case, 
the molecule) , defined as the difference between the total 
entropy S(T) and the entropy calculated in the absence 
of the molecule S*- ' (T) . As it has been shown by Bcthc- 
ansata^ and NRG calculations^ in two-channel Kondo 
models, S'imp (T) = S(T) - S^{T) reaches a universal 
low-temperature plateau at Si mp /kB = l/21n(2) {ks is 
Boltzmann's constant) when both channels are equally 
coupled to the impurity. 

Fig. [5] shows Si m p versus temperature T/D and a for 
two values of the molecular level energy: Sd = —U/2, cor- 
responding to the particle-hole symmetric point (top pan- 
els), and Ed = — 0.9J7 (bottom panels).— In both cases, 
an entropy plateau Simp/fcs = 1/2 ln(2) is reached at low 
temperatures (Figs.[2ja) and (b)), signaling the presence 
of a non- Fermi liquid (NFL) fixed point. 

This NFL fixed point is reached as a approaches a 
critical value a c and at temperatures below a character- 
istic crossover energy scale T*. Deviations from a — a c 



FIG. 3: (color online) Critical dot-lead coupling a c (a) and 
crossover temperature T* (b) vs. molecule level position Ed- 
(c,d) Same quantities but now as a function of the electron- 
phonon coupling A at Ed = —U/2. 



drive the system away from this state to the more conven- 
tional Kondo-screened state (characterized by S'imp = 0) , 
illustrating the unstable nature of the fixed point. This 
is depicted in Fig. H{c), showing S; mp versus a at fixed 
temperatures. At lower temperatures, a narrow peak of 
height S; mp = 1/2 In 2 pinpoints the critical value a = a c . 
This marks the position of the fixed point in parameter 
space. At higher temperatures, the broadening of these 
peaks indicate that signatures of NFL behavior extend 
over a finite range of a. 

This behavior persists away from particle-hole symme- 
try, as shown in Figs. HJb) and (d). Interestingly, the 
critical value a c increases as the system approaches the 
mixed-valence regime (e<j — * 0, —U ; (rid) — ► 0, 2) as com- 
pared with the Ed = —U/2 case {{rid) = !)■ This con- 
trasts with the (phonon-independent) two-channel An- 
derson model (2chAM), for which the critical couplings 
are erf-indcpcndenti^ The crossover energy scale to the 
NFL fixed point, related to the Kondo temperature, de- 
creases exponentially as the system enters the mixed- 
valence regime, as in the 2chAM. 

In the mixed-valence regime, the system flows into a 
Fermi-liquid fixed point with Si mp = fcs hi(2) at higher 
temperatures before entering the NFL regime (as seen, 
e.g., in Figf^b)). This fixed point is nonmagnetic and 
characterized by a parity degeneracy in the ground state 
(rather than the usual spin degeneracy in similar models). 

As previously discussed, the Si mp vs. T/D curves dis- 
play the NFL plateau at low temperatures for — U < 
Ed < at a critical a = a c . The critical value a c in- 
creases as the system is moved away from particle-hole 
symmetry, reaching its highest values at Ed — 0~, — U + , 
as depicted in Fig. [3Ja). In addition, at a — a c the 
crossover temperature sharply decreases as the systems is 
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FIG. 4: (color online) (a) Conductance G versus Ed for dif- 
ferent values of a. A dip appears for a = a c at Ed — — U/2, 
signaling the 2chK state, (b-d) Even (S e ), odd (S ) and rela- 
tive \8 e — 8 \ phase shifts vs. temperature at Ed = — U/2 and 
different values of a. 



driven away from the particle- hole symmetric point. This 
is illustrated in Fig. [3jb) , where we define the crossover 
temperature T* as S imp (T = T*) = 3/4 In 2. This panel 
shows that T* decreases exponentially as Ed —* —U, 0, in- 
dicating that the 2chK is only reached within the range of 
gate voltages for which the molecular level in singly occu- 
pied. In fact, we find the NFL fixed point in this "local 
moment" range (-U < Ed < 0) only, indicating, along 
with susceptibility calculations, that the Kondo screen- 
ing occurs in the "real spin" as opposed to a "pseudospin" 
degree of freedom (e.g., those connected to charge states 
or \m),\m ± 1) phonon states). In this range, the value 
a c varies nearly quadratic with Sd, as shown in Fig. |5fa). 

Furthermore, for fixed £<j, both a c and T* decrease 
as the electron-phonon coupling A increases (Figs. [3]-(c) 
and (d)). The decrease in T* is consistent with the 
fact that the electron-phonon coupling renormalizes the 
electron-electron interaction as U c g « U — 2A 2 /wo^ effec- 
tively driving the system into the mixed valence regime. 
The decrease in a c with A indicates that not only the 
molecule-lead couplings but also the coupling of the vi- 
brational mode with the charge state of the molecule 
plays a role in the mechanism leading to the two-channel 
Kondo effect. 



IV. PHASE SHIFTS AND CONDUCTANCE 

We now turn to the transport properties across the 
molecule. We first point out that the presence of the 
CM phonon term in Eq. (JTJ breaks the "proportionate 
coupling" condition^ between left and right leads (even 
in a rather unrealistic symmetric coupling configuration 



Vl = Vr). For this reason, a calculation of the linear con- 
ductance via a Landauer-like formula 2 ^ would involve not 
only the local interacting retarded Green's function (ob- 
tainable from NRG) but also Keldysh Green's functions. 

Instead, we turn to the equivalent approach of calcu- 
lating the T — conductance using the scattering phase 
shiftS ) 14 i 21 ' 27 as G = Go sin 2 (5 e — 5 a ), where S e and S 
are the phase shifts in the even and odd channels, re- 
spectively, and Go = (2e 2 /h) sin 2 9 with 9 being an over- 
all phase that depends on the microscopic details of the 
molecule- lead junction (we henceforth consider 9 = n/2). 

The phase shifts S e t a \ can be obtained from the NRG 
spectr a 28 i 29 using the (Q, S, P) quantum numbers to la- 
bel the states^ We note that the parity quantum num- 
ber is strictly conserved only for A = 0. Although it is 
possible, in some cases, to calculate the difference 5 e — 8 
from the (Q, S) NRG spectra ) 14 i 29 in the following we use 
A = as it retains most of the interesting physics . We 
should note that, away from FL fixed points, the corre- 
spondence between the excitations in the NRG spectra 
and the phase shifts entering the conductance formula 
is only approximate. Nevertheless, we expect the con- 
ductance obtained with this prescription to give a qual- 
itatively accurate picture in the NFL state as well, as 
discussed below. 

Results for the conductance are shown in Fig. HJa). 
For a = 0, the familiar shape is recovered: G = Go at 
the particle-hole symmetric point Ed = — U/2 an G — ► 
as \ed\ increases. As a increases, the peak narrows with 
G = Go at the p-h symmetric point for a < a c and 
a > a c . Interestingly, as a approaches the critical value 
a c at Sd = —U/2, a dip appears in the conductance curve. 

This is a indication of the NFL behavior and a signa- 
ture of the two-channel fixed point. The behavior of indi- 
vidual phase shifts at the particle-hole symmetric point 
the is illustrated in Figs(3Ib,c,d). For a < a c (FigHJL)), 
5 e — ^ 7r/2 and 5 a — > at low temperatures, indicating a 
decoupling of the odd channel, while for a > a c the odd 
channel becomes strongly coupled at low temperatures 
(S — > 7r/2, S e — » in Figs|3]-(d)). In both cases, one 
expects a peak in the conductance^ 

At the critical point a — a c (NFL regime), the NRG 
spectra is identical for both even and odd parities, as 
predicted by conformal field theory^ In this case, our 
prescription for obtaining the phase shifts gives S e = 5 Q 
at low temperatures, as depicted in Fig. Hfb), causing 
a destructive interference and suppressing the transmis- 
sion. We thus expect this result to hold, even though the 
individual values of <5 e ( D ) obtained from the NRG spectra 
in the NFL regime are only approximate. 

These signatures in the low temperature conductance 
versus gate voltage curves can, in principle, identify the 
2chK regime in molecular junctions. Experimentally, a 
fine tuning of the microscopic parameter a to the critical 
value is nonetheless a challenging task. In general, the 
value of a will be determined by specific details of the 
junction, such as the ratio of amplitude of the center- 
of-mass vibration and the length of the molcculc-lcad 
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overlap. It is interesting, however, to notice that, at 
higher temperatures, one might obtain S e ~ S over a 
wider range of a and Sd near the critical values, as il- 
lustrated in Figs. [2)Jc)-(d). In this case, one expects a 
finite-temperature signature of the T = non-Fermi liq- 
uid point. 

V. SUMMARY 

In summary, we have studied center-of-mass vibra- 
tional effects and phonon-assisted processes in the trans- 
port properties of a molecular junction in the Kondo 
regime. The interplay between electron-electron and 
electron-phonon interactions in this system can be de- 
scribed by an effective two-channel Anderson model 
with phonon-assisted tunnel couplings. Our numeri- 
cal renormalization-group calculations for the thermody- 
namic properties of the effective model show non-Fermi- 



liquid effects below a characteristic crossover tempera- 
ture over critical lines in parameter space. 

We find that the crossover temperature is at a max- 
imum at the particle-hole symmetric point and rapidly 
approaches zero as the system enters the mixed valence 
regime. Furthermore, we find distinct signatures of the 
non-Fermi-liquid phase in the linear conductance. 
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